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Abstract 

For cosmic neutrinos we study the conditions and the effects of the 
coherence loss as well as coherent broadening of the spectrum. We 
evaluate the width of the neutrino wavepacket produced by charged 
particles under various circumstances: in an interaction-free environ- 
ment, in a radiation-dominated medium (typical of the sources of the 
gamma ray bursts) and in the presence of a magnetic field. The ef- 
fect of the magnetic field on the wavepacket size appears to be more 
important than the scattering. If the magnetic field at the source is 
larger than ~10 Gauss, the coherence of neutrinos will be lost while 
traveling over cosmological distances. Various applications of these 
results have been considered. We find that for large magnetic fields 
(B > 10 9 Gauss) and high energies (E u > PeV), "coherent broad- 
ening" can modify the energy spectrum of neutrinos. In the coher- 
ent case, averaging out the oscillatory terms of the probabilities does 
not induce any statistical uncertainty beyond what expected in the 
absence of these terms. A deviation from the standard quantum me- 
chanics that preserves average energy and unitarity cannot alter the 
picture. 

PACS: 14.60.Pq 

KEYWORDS: Cosmic neutrinos, Coherence, Spectrum 



1 yasaman® theory, ipm . ac . ir 

2 smirnov@ictp . it 



1 



1 Introduction 



Construction and operation of the neutrino telescopes open new windows 
towards exploring the Universe. By studying the neutrinos arriving at these 
telescopes from sources located at cosmological distances, we can derive un- 
precedented information on the sources of such neutrinos as well as on the 
propagation and properties of the neutrinos, themselves. 

Operating detectors have not so far detected cosmic neutrinos. The data 
accumulated during 2000-2004 by the AMANDA experiment shows no indi- 
cation of neutrino flux from beyond the atmosphere of the Earth [1] which 
leads to an upper bound on the point source flux 

£ 10- 8 GeV cm- 2 sec" 1 sr" 1 , (1) 

dE u 

in the energy range between 1.6 TeV and 2.5 PeV. This bound implies that 
the total number of events per year at a 1 km 3 scale detector such as ICE- 
CUBE [2] cannot be on average larger than 10 4 . A host of theoretical models 
predict cosmic neutrino fluxes in the (1 -100) TeV energy range which can be 
detected by a 1 km 3 scale detector: e.g., the fireball models for Gamma Ray 
Burst (GRB) sources [3, 4], models for supernovae type I b/c [5, 6] and mod- 
els for the Active Galactic Nuclei (AGN) [7]. The neutrino flux from GRB 
in the (1-10) TeV range can saturate the present bound from AMANDA [4]. 
We should prepare for such generosity of nature and anticipate what can be 
learned about the neutrino properties as well as the mechanism behind the 
production of neutrinos at the source. Because of this prospect, we mainly 
concentrate on neutrinos with energies in the (1 - 100) TeV range throughout 
the present paper. 

According to the models, neutrinos are produced through the following 
chain of reactions [8] 

p + 7 -> A -> 7T+ + X , (2) 

fi + -> v e + + e + . (3) 

Along with the process (2), both n + and tx~ can be produced via the pp colli- 
sions or via the collision of the protons with the nuclei present in the medium. 
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If the interaction rates of tt and [i are negligible, the flavor composition of 
the neutrino flux at the source will be as the following 

+ = (K + F °») ■ + = 1 : 2 : , (4) 

where F° and respectively denote the fluxes of v a and V a . 

Neutrino oscillations alter the flavor composition of the neutrino flux dur- 
ing propagation from the source to the detector. If the flavor composition 
at the source is known, by studying the flavor composition at the detector, 
one can extract information on the oscillation parameters [9]. In particular 
in [10, 11], as the key part of the program of reconstructing the unitarity 
triangle and deriving the lepton sector Jarlskog invariant, it was suggested 
to study the cosmic neutrinos to extract the U^i element of the Pontecorvo- 
Maki-Nakagawa-Sakata matrix. There is also rich literature on the possibility 
of directly determining the Dirac CP-violating phase from the neutrino tele- 
scope data [12]. 

One of the important effects of neutrino propagation which can in prin- 
ciple influence observations is the loss of coherence. The different neutrino 
mass eigenstates having the same energy have different velocities. Thus, the 
wavepackets of the mass eigenstates composing a neutrino state will come 
apart as they propagate. If the traveled distance is so large that these com- 
ponents completely separate from each other, they will cease to interfere at 
the detector. In the case of coherence loss, the oscillatory terms of the oscil- 
lation probabilities disappear. Conversely, in the coherent case, the different 
mass eigenstates keep interfering and as a result, the oscillatory terms are 
still present. Of course, even in the latter case, because of the finite energy 
resolution of the detector and the fact that the sources of different events are 
located at various distances, the effects of the oscillatory terms average to 
zero. As discussed in [13, 14], the two cases are in practice indistinguishable. 
In this paper, we discuss in detail under which conditions the coherence is 
lost. We also evaluate the statistical uncertainty in the presence of the oscil- 
latory terms for the coherent case. 

Coherent broadening of the wavepacket can lead to deformation of the 
energy spectrum which has direct observational consequences. In this paper, 
we discuss under what circumstances, the effect will be noticeable. 

The paper is organized as follows. In sec. 2, we elaborate on the conditions 
of the coherence loss. In sec. 3, we explore how the interactions of the 
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parent particles in the production region can affect the sizes of the emitted 
wavepackets, and consequently, the loss of coherence. In sees. 4.2 and 4.3, 
we discuss the applications of our finding for neutrinos from various cosmic 
and terrestrial sources. Sec. 4.4 is devoted to the discussion of coherent 
deformation of the energy spectrum. In sec. 4.5, we consider the effects of 
incomplete averaging. In sec. 4.6, we comment on how a deviation from the 
standard quantum mechanical evolution of states can affect the results. A 
summary and the discussion of results are given in sec. 5. 



2 Loss and restoration of coherence 
2.1 Coherence loss 

The oscillatory terms from the oscillation probabilities can disappear in two 
physically distinct situations: (1) loss of coherence; (2) averaging. The first 
case is related to features of production of neutrinos, whereas the second one 
- of detection. For the scope of further discussion we first formulate the 
condition for loss of coherence in the momentum space and then derive the 
same condition considering the problem in the configuration space. The two 
approaches are equivalent and give the same result. 
The phase of oscillation, 0, equals 

«I5,L) S ^, (5) 

where Am 2 is the mass squared difference and L is the distance between the 
source and the detector which hereafter will be referred to as the baseline. Let 
us define AE L as energy difference for which (f>(E — AE L , L) — <f)(E, L) = 2n. 
It is straightforward to show that 

where l v = AttE/ Am 2 is the vacuum oscillation length. Neutrinos traveling 
over cosmological distances are essentially on their mass shells. We can there- 
fore use the dispersion relation: E 2 = p 2 + m 2 . Introducing the uncertainty 
in the energy, oe and in the momentum 3 , a p , we find from this relation 



3 From now on, by a p we implicitly mean the width of the wavepacket in the longitudinal 
direction. 
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that Eoe = p<J p . Width of the wavepacket in the energy-momentum space is 
determined by <je and a v . Furthermore, since neutrinos are ultrarelativistic 
(E — p), the uncertainty in energy and momentum are practically equal: 
a E ~ Up. In what follows we will not distinguish between a E and a p . 4 

The interference between the effects of different mass eigenstates disap- 
pears if 

«» ^ = = 4 (7 » 

(For an explicit proof of the disappearance of the oscillatory terms see the 
appendix.) Omitting the numerical factor, 4-7T, we can rewrite the condition 
in Eq. (7) as 

Am 2 L a v _ . . 

— E >y (8) 
Let us now reconsider the same problem in the configuration space. The 
size of the wavepacket in the configuration space at the production point is 
given by the inverse of the uncertainty a p : 

a x ~ — . (9) 

A neutrino of a definite flavor can be considered as a superposition of the 
wavepackets of different mass eigenstates. Due to the difference between the 
group velocities of the mass eigenstates, 

Am 2 . . 

= w (10) 

the corresponding wavepackets separate in the course of propagation. The 
separation, d^, after propagating a distance of L is equal to 

L . Am 2 

The separation in the configuration space is complete when 

An~i 

a x « d L = L—. (12) 



4 Under certain circumstances, like emission of neutrinos from crystals with Mossbauer 
effect, the difference between the energy and momentum uncertainties becomes important. 
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At this point, the overlap of packets and therefore interference effects dis- 
appear - the coherence is lost. The condition in (12) is equivalent to the 
condition in (8) in the energy-momentum space. The inequality (8) as well 
as (12) is called the condition for "loss of coherence". The coherence length 
is defined as 

E 2 

L coh = 4 ™*^- (13) 

For L > L coh , the coherence is lost and if the coherence is not restored at 
the detection (see subsec. 2.2 as well as the appendix), the oscillatory terms 
disappear from the oscillation probabilities. 

For the typical source to detector distance L = 100 Mpc, the neutrino 
energy E — 10 TeV and Am 2 = Am 2 atm = 2.5 • 10~ 3 eV 2 , from Eq. (12) we 
find that the coherence is lost if the size of the wavepacket satisfies 

, o ,n a ( L \f Am2 WlOTeV\ 2 
a x < d L = 3 • 10~ 3 cm — - — ? — . 14 

VlOOMpcy V2.5-io- 3 ev 2 y V E ) K J 

Replacing Am 2 with Am 2 oZ , we obtain a value for di that is about 30 times 
smaller. That is, if the size of the wavepacket is much smaller than ~ 10~ 4 
cm, the oscillatory terms given by Am 2 ol disappear. 

In our consideration we have ignored spread (widening) of each wavepacket 
in the configuration space due to the presence of different energies/momenta 
in the packet. For a given mass eigenstate, the group velocities of the compo- 
nents of the wavepacket with AE = <je are different: Av = {m 2 / E 2 ){a e / E) . 
The spread of the wavepacket over a distance L is therefore equal to 

ffL= |A.L| = f!fi=e) 3 ^L. (15) 
E 2 E \EJ a x m 

Comparing this expression with the separation (11) we find that the condition 
d-L °l implies 

o E Am 2 

-T7 « — 

E m A 

where m is the heaviest neutrino mass. The minimal value on the RHS 
corresponds to the degenerate neutrinos. Eq. (16) is fulfilled for all situations 
under consideration in this paper. That is, the separation of the wavepackets 
is more important than their spread. 



4 « — r, (16) 
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2.2 Restoration of coherence 

Separation of the wavepackets does not mean yet that the effects of oscil- 
latory terms are absent in observation. If the time scale of the detection 
process is larger than the time interval between the arrival of the successive 
wavepackets composing a single state, the wavepackets still interfere at the 
detector despite the spatial separation [13]. In other words, if the detector 
has large enough "time memory" the detection restores the coherence. The 
coherence is determined both by the production and detection. In this sense, 
there is a symmetry between the source and detector. 

Consider a neutrino wavepacket whose components are completely sep- 
arated: i.e., a x ~ dp 1 <C d L , where d^ is given in Eq. (11). If the en- 
ergy of the neutrino is measured with precision better than dj^ 1 , the coher- 
ence will be restored [13]. According to the energy-time uncertainty princi- 
ple, such a measurement takes a time longer than d^jc during which the 
second component arrives at the detection point. However, in the re- 
alistic situations under consideration, restoration of the coherence at de- 
tection, which requires determination of energy with precision better than 
d~ L l ~ 1(T 2 eV(100 Mpc/L)(2.5xl(r 3 eV 2 /Am 2 )(£/10 TeV) 2 (see Eq. (11)), 
does not take place. Suppose that a hypothetic detector reaches such a sen- 
sitivity; then, for each individual neutrino state, coherence is restored. How- 
ever, even in this, once we accumulate data from different sources located at 
various distances, the oscillation pattern will disappear. 

3 Width of the neutrino wavepacket 

In this section, we evaluate the width of a neutrino wavepacket produced 
through the processes shown in Eq. (3) in an environment such as the rela- 
tivistic jets which, according to the fireball models, are present in the GRB 
sources. We discuss the dependence of the wavepacket size on the density 
and the strength of the magnetic field at the production region. We also 
discuss if a neutrino produced in such an environment loses its coherence 
traveling over cosmological distances. 
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3.1 Neutrino produced by decay of "free" pions and 
muons 

In this subsection, we consider the case in which the parent charged particles 
(pions and muons) do not undergo any interaction with matter or with the 
magnetic field before decay. The consideration applies for any environment 
in which the mean free path and the Larmor radius of the parent particle are 
much larger than the distance that the charged particle travels before decay. 
The cosmogenic neutrinos (i.e., the neutrinos that are produced through the 
interaction of the cosmic rays with the galactic and intergalactic medium) 
can be considered as an example for realization of such a situation. 

The width of the wavepacket of the neutrino in the rest frame of the 
parent particle can be estimated as 

o° E « a° p ~ - , (17) 

TO 

where r is the lifetime of the parent particle in its rest frame (i.e., for the 
muon r = 2.2 x 10~ 6 sec while for the pion r = 2.6 x 1CT 8 sec). 

Suppose the pion propagates along the z-direction; i.e., in the observer 
frame, the energy-momentum of the pion can be written as (E n , 0, 0, p n ). Let 
us take the four-momentum of the neutrino in the rest frame of the pion to 
be 

(£°,0,£°sin£,£°cosO, 

where 

„ ml — ml , 

^•-kr-"-' 4 (18) 

and £ is the angle between neutrino momentum in the frame of pion and 
the momentum of pion in the observer frame. The energy of the neutrino 
in the observer frame can be then written as E v = (E n /m w )E®(l + v n cos£) 
where v n = p^/E^ ~ 1 is the velocity of the pion. Performing a Lorentz 
transformation we find that in this frame 

E E AE 1 
oe = — (<J E + v n a cos£ = • 19 

Notice that oeIE v is Lorentz invariant. 

It is instructive to derive these results through alternative approaches. 
The length of the wavepacket in the rest frame of the pion is equal to cr . It 
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is straightforward to show that in the observer frame the wavepacket length 
is Lorentz contracted to 

7(1 + cos 4) E u AE U 

The last equation is for neutrinos emitted in the forward direction £ = 0. As 
expected, Eqs. (19,20) yield a x = cr^ 1 . 

We can rederive Eq. (20) by using the fact that the number of wave 
periods in the wavepacket is a Lorentz invariant: 

— = o x E v = invariant. 
A 

In the rest frame of the pion this number is r°p° = t°E®, so 

ax = ~ r ^. (21) 

E v AE V 1 ; 

Neutrinos emitted in the forward direction have higher frequency (energy) 
but shorter packet (pulse). 

Another Lorentz invariant that one can use to determine a x in the ob- 
server frame is the number of the oscillation periods before the system loses 
its coherence. This number is given by the ratio of the coherence length L co h 
to the oscillation length, l v . Using Eq. (13), we find 

— ^ = a x E v = invariant, (22) 

which immediately leads to the same result as (21). 

Numerically, for neutrinos produced by pions we obtain 

a x ~ 2 x 10~ 3 cm ^ 10 ^ eV ^ . (23) 

For the neutrinos produced by muons, o x is larger by the ratio of the lifetime 
of the muon to that of the pion: a x ~ 0.2 cm(10 TeV/E u ). 

Let us consider the ratio of separation length and the length of the 
wavepacket which, as discussed in sec. 2.1, determines the loss of coherence: 

d L _ Am 2 L 1 Am 2 L 2 

V x ~ 2ElV x ~ E v m^' { } 
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(For neutrinos from the muon decay one should substitute m n by and 
insert the lifetime of muon for To.) Notice that this ratio depends on the first 
power of the neutrino energy. Numerically, we find 

<'l _ , ( Am' N f L \ /lOTeV\ ( 3xUT^ec\ 

~ x ~ 01 UxlO-eVj [woWc) \-^r) { 7 J ' (25) 

This ratio shows that despite the extremely long baselines, coherence is in 
general maintained. This is always correct for neutrinos with E v > 1 TeV 
originating from the muon decay (even if we take L as large as the present 
size of the Universe). Only for relatively low energy neutrinos, E v ~ TeV, 
from the pion decay, the coherence can be lost. 

3.2 Interaction with particles in medium 

If the parent particle (pion, muon) scatters off the particles in the medium, 
the length of the wavepacket of the produced neutrino will be shorter than 
in the interaction-free case. Here, for illustrative purposes we consider a 
radiation dominated environment as expected in most models for the GRB 
sources. According to the fireball models [4, 15], neutrinos at the GRB 
sources are produced in relativistic jets which have boost factors (Tj et ) of 
order of 10 2 in the observer frame. In the frame comoving with the jet, 
which we will call the jet-frame, photons have a thermal distribution with 
a typical temperature of T 7 ~ few keV [4]. In this environment, the main 
interaction process for the charged particle is scattering off the photons. 

Let us consider the relevance of collisions and obtain the size of the neu- 
trino wavepacket a J x et in the jet-frame. The most energetic pions (and muons) 
are supposed to move in the direction of jet and therefore their energy in the 
jet frame can be estimated as 

^^ 100G ev(^)(g), <*> 

where E n {E^) is the energy of the pion (the muon) in the observer frame. 

The average distance that the parent particle travels between two succes- 
sive collisions can be estimated as 

4oi ~ — , (27) 
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where n 7 = [2£(3)/-7r 2 ]T 3 is the photon density and a is the total cross-section 
of the scattering. In the photon-pion center of mass system, the momentum 
of the pion is of order of ~ ElT^/m^ ~ 4 MeV <C m^; that is, in this frame 
the pion is non-relativistic. In this regime, the cross-section of the Compton 
scattering is given by the Thompson formula: 



In general, to calculate £ co i, we should take into account the energy distri- 
bution of the photons in the thermal bath but since in this non-relativistic 
regime the cross-section is not sensitive to the momentum of the photon, the 
approximation we have made is justified. Inserting n 7 and a in Eq. (27), we 
obtain 



In the case of muon scattering, m n has to be replaced by m M . Notice that 
the mean free path is approximately the same for muons and pions. 

In the jet-frame the pions are ultra-relativistic so the majority of the 
produced neutrinos are emitted within a narrow cone with an opening angle 
of m n / El, oriented along the momentum of the parent particle. We will refer 
to this cone as "the emission cone." The energy transfer in each collision is 
negligible (i.e., the collisions are elastic) and only the momentum of the 
charged particle randomly changes direction by an angle of size ~ T^/m n . 
The smallness of the momentum change in each collision (i.e., T y /m n <C 
m n /El) implies that if the direction of the momentum of the parent particle 
before collision is such that the line of sight lies within the emission cone, 
after the collision, the line of sight will still remain within the emission cone. 

All the consideration we had so far in this sub-section holds equally for 
the pion and the muon. However, we should bear in mind that the neutrino 
production by the former takes place via a two-body decay while from the 
latter it takes place via a three-body decay. As a result, the widths of the 
neutrino wavepackets produced by them are different. In the following we 
first discuss the pion decay and we then consider the muon decay. 

During the time between two successive collisions (i co i/c), the pion emits 
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a X,lT ~ Zcol — — = Z-COI — — — , (29) 



a wavepacket whose size is of order of 

m^-Eg = m n ml - 
EiEl c ° l El 2m n Ei 

where E® is the energy of the neutrino in the rest frame of the pion (see 
Eq. (18)). Notice that, in Eq. (29), we have taken into account the Lorentz 
contraction of the wavepacket as discussed in sec. 3.1 (see Eq. (20)). 

After the collision, the momentum of the pion is changed by \AfP w \ ~ 
TyE^/m n in a direction transverse to the momentum (Ap£ - pi -C |Ap£| 
Writing the kinematics of scattering, we find that the difference between the 
energy of the neutrino emitted in our direction before and after the collision 
is ~ T^Ei/m^) 2 which is much larger than (cr^*)" 1 given in Eq. (29). This 
means that the wavepackets emitted before and after the collision are inco- 
herent and the effective size of the wavepacket in the jet frame is therefore 
given by Eq. (29). 

The spectrum of the neutrinos produced by the muon is continuous. Thus, 
even after a collision the wavepacket of the neutrino will have a component 
of the same energy as before. However, if we wait long enough the direction 
of the parent particle will change so much that the line of sight will exit 
the emission cone. After this time, the amplitude of the neutrino emitted in 
the observation direction is negligible. Thus, the size of the wavepacket is 
determined by the time, At, during which the line of sight exits the emission 
cone. Let us estimate At. As mentioned, after each collision the momen- 
tum of the parent particle changes direction randomly by an angle of size 
T 7 /m M so, after N collisions, the total rotation angle on average amounts to 



NT 7 /m^. Thus, for E^ < mj, after N ~ [m^jE^T^ 2 collisions, the 
direction of the momentum changes by an angle of ~ m^/E^ and the line 
of the sight therefore exits the emission cone. Since the time between two 
successive collisions is of order of £ co i / c, we conclude that 

2 



At 



<c„l / "if, 



^T 7 



Taking into account the Lorentz contraction (see sec. 3.1) for neutrinos emit- 
ted in the forward direction we obtain 

<~cAt^g~% < , (30) 
^ El El 3 T 2 (El) 2 E 2 ' y 1 
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where E° is the energy of the neutrino in the rest frame of the muon which 
is of order of m M /3. Notice that this relation is valid only if E^T^ <C m 2 . For 
larger values of EjjTj, a single collision will be enough for the line of sight to 
exit the emission cone (i.e., N = 1) and therefore ~ £ co i(m^E°) / (E^El) . 

As we discussed in the previous section, <je/E, or equivalently Ea x , is 
invariant under the Lorentz transformations. Thus, from Eq. (29), we find 
that in the observer frame, the wavepackets of the neutrinos produced by the 
pion are of size of 

With this wavepacket size, we obtain from Eq. (14) that for neutrinos pro- 
duced by the pions, the wavepackets of all three neutrino mass eigenstates 
will be separated on their way to Earth from the GRB sources. 
Similarly, using the invariance of Ea x we obtain 

„ mt = a ^ „ 5 x 10- c m (i™)' (i°™V (Em) 3 (32 ) 

* E v V T 7 / V Ev ) V 100 J K ' 

From Eq. (14) we conclude that unlike the case of neutrinos from the pion 
decay, the different mass components of the neutrinos produced by the muon 
will keep overlapping until they reach the Earth. In other words, the coher- 
ence of the neutrinos from the muon decay is maintained. 

In the energy-momentum space, the wavepacket sizes in Eqs. (31,32) 
translate into 

and 

256 „ 10- (-*L-) '(-$-)' (™) ' (£or Ml « m l). (34) 

Increasing the energy, oe increases. Notice however that for E V T~,/T ' Je t ~ m 2 , 
the above formulas are not valid for two reasons: (i) in this regime the cross 
section is not given by the Thompson formula; (ii) as discussed before, for 
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high energies, Eq. (30) is not valid. Taking into account these considerations 
we find that 

~^~5x 1(T 12 (-^-) 3 (for ^ £ ml). (35) 



E„ E v MkeV v T 



jet 



3.3 The effects of magnetic fields 

Let us discuss the effects of the magnetic field present in the source on the 
width of the neutrino wavepacket. The magnetic field at the GRB sources 
can be as large as ~ 10 9 Gauss [4, 5, 6]. For simplicity, we assume that the 
magnetic field is constant and uniform: (B = Bz). Later on, we will show 
that this assumption is justified. Notice that although the magnetic fields in 
the sources under consideration are quite large but still eB <C E 2 and we can 
therefore treat the effects of the magnetic field classically. That is, we can 
describe the state of the charged particle as a plane wave whose momentum 
in the direction transverse to the magnetic field slowly rotates in time. As 
we will see, the decays of the pion and muon have to be treated differently. 

Let us first discuss the pion decay. In the magnetic field, the charged 
pions move on circular (spiral) trajectories with Larmor radius R = E n /eB. 
The energy of the emitted neutrino in a two body decay depends on the angle 
between its momentum and the momentum of the pion which, in turn, is a 
function of time. That is, the energy of the emitted neutrino at time t, E u (t), 
differs from that at time t + At: 

AE v = E v (t + At)-E u (t). (36) 

Taking into account the Lorentz contraction (20) we find that the length of 
the neutrino wavepacket in the configuration space emitted during the period 
At from a pion with energy of E n is given by 

a x = (cAt) — — , 

where E® is the energy of the neutrino in the rest frame of the pion (see 
Eq. (18)). Consider segments of the pion trajectory for which AE V defined 
in Eq. (36) coincides with oe = v^ 1 : 

GE = ak = ai © (§) = it (^^|) • (37) 
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The wavepackets emitted during successive At can be resolved by a hypo- 
thetical detector with energy resolution <je and time resolution a^ 1 , so their 
effects sum up incoherently. Consequently, the coherence length of the neu- 
trino wavepacket is determined by cr^ 1 which can, in turn, be derived by 
solving Eq. (37). For this, one needs to determine AE U or E v as a function 
of time. 

The four-momentum of the pion can be written as 

(E n , p n sin 6^ cos p^ sin 6 n sin &(t), p n cos9 n ), (38) 

where <3>(i) is the rotation phase: 

= S^ZIA . (39 ) 

Notice that E n is constant. Let us denote the unit vector in the direction 
from the source to the detector by /, and define the axis y as y — I x z/\l x z\. 
In this frame, the neutrino emitted at time t towards the detector has the 
following four-momentum 

(1, sin 0„, 0, cos (40) 

where cos V — (l • z) and 

EM = '"- " < (41) 

2[E n — p-x cos(6 u — 9 n ) + 2pn sin 9 n sin 9 V sin 2 -^] 

Notice that due to the two-body character of the pion decay, the neutrino 
energy in the rest frame of the pion is fixed, and in the observer frame it is 
uniquely determined by the angle between the momenta of the pion and the 
neutrino. From Eq. (41), we obtain 

AE U = 2e^sm^sin^sin$ At _ 
m\ — m^ 

Only a small fraction of the neutrinos (of order of (m^/p^) 2 ) are emitted 
in directions for which the angle between the neutrino and pion momenta 
is larger than m^/p^. We therefore concentrate on the angles and emission 
time for which the line of sight lies inside the emission cone: 

(0„ - 6 V ) ~ sin $(t) ~ mjp v < 1. (43) 
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For such values of t, 9 n and V , from Eqs. (37,42) we obtain 

4eBE^m n sin 9 n sin 6 U 1 1 ' " 



(ml - m 2 ^) 2 



(44) 



Let us apply these results to neutrinos which are produced in relativistic 
jets [4, 15]. We perform our analysis in the (comoving) jet-frame where the 
electromagnetic field is predominantly composed of the magnetic component 
5 . Using the invariance of &e/E (see Eq. (19)), we can immediately write 
the width in the observer frame, oe, in terms of the energy width in the 
jet-frame, o 3 ^: 



oe = o^_ / eBE„m n 4 / 100 E u B \ 1/2 

E v ~ y T jet (ml - ml) 2 ~ \T iet 10 TeV 10 7 Gauss J ' 

(45) 

where we have taken 4 sin 6 n sin 6 U ~ 1 and E ] v et ~ E u /Tj et . The width of the 
wavepacket in the configuration space is given by 

-i o in-" f r jet W 7 Gauss V 7 V 10 TeV V /2 
^~^e ~2x 10 cm -j- — . (46) 



100 B ) \ E v ) 

Conferring the estimations in Eqs. (46) and (31) we conclude that for B > 
10~ n Gauss, interaction of pions with the magnetic field is more effective 
in shortening the wavepacket than their scattering off the particles in the 
medium. That is, in practice, the size of the wavepacket of neutrinos pro- 
duced in a cosmic neutrino source is given by Eq. (46) rather than by Eq. (31). 

Using the ratio in Eq. (45) and the condition shown in Eq. (8), we find 
that the loss of coherence over cosmological distances occurs if 

.„ „ f8 x 10~ 5 eV 2 \ Vl00 Mpc\ 2 E v Y jet 



B»5x 10- 13 Gauss \ 9 r ^ (47) 

V Am 2 J V L J 10 T eV 100 v ; 

Let us now consider neutrinos from the muon decay. Due to the three- 
body character of the decay, the energy spectrum of the neutrino emitted 
in a given direction is continuous. Thus, the consideration that we had in 



5 In the observer frame, according to the Lorentz transformation in addition to the 
magnetic field, an electric field will be also present. 
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the case of the two-body decay of the pion does not apply here. Since the 
muon is ultrarelativistic, the majority of the emitted neutrinos are oriented 
in the direction of the muon momentum within the emission cone. Only a 
small fraction of 0{m L jE' l \ is emitted outside the emission cone. Since the 
momentum of the muon rapidly rotates, only for a short period, At, the line 
of sight lies inside the emission cone. A reasonable estimate for At is the time 
interval during which the angle between the line of sight and the momentum 
of the muon is smaller than m^/E^. (eBAt/E^) ~ m^/E^. We then find 
that in the jet-frame 

^(aW<V = ^M 2 , (48) 



E 



where we have taken E u ~ E^/3. Using the invariance of (Je/E, we obtain 

= Ejt = 5 x lO" 8 — ^ ^—^ (49) 

E v Ei et 10 7 Gauss 10 TeV T jet ' 1 ; 

where E v is the energy in the observer frame: E v = Ei et /Tj et . Thus, in the 
observer frame, the length of the wavepacket is 

i K in -n f 10 TeV V 10? Gauss Tjet 

Again, comparing the estimations in Eq. (50) and in Eq. (32) we find that 
for B > 10 _1 Gauss, the interaction of muons with the magnetic field is 
more effective in shortening the wavepacket size than their scattering off the 
particles in the medium. In the presence of the magnetic field, the wavepacket 
size is therefore given by Eq. (50) rather than by Eq. (32). 

To obtain complete loss of coherence (see, Eq. (8)), the following condition 
has to be fulfilled: 

100 MpcSxlO- 5 eV 2 T jet 
B » 10 Gauss— ^ lOO"- (51) 

Notice that while in the case of the pion decay, a E is proportional to B 1 ^ 2 
(see Eq. (45)), in the muon decay case, the dependence on B is linear (see 
Eq. (49)). This explains the huge difference between the lower limits on B 
shown in Eqs. (47,51). 
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In the analysis, we have taken the magnetic field to be constant. In fact, 
this assumption is justified as long as the variation of the magnetic field 
over scales of cAt [where At is the time segment that determines <7e\ see 
Eqs. (37,48)] is negligible. Taking E v = 10 TeV and T jet = 100, we find that 
for the case of pions, cAt ~ 10~ 5 cm(10 7 Gauss/S) 1 ^ 2 and for the case of 
muon, cAt ~ 4 x 10~ 3 cm(10 7 Gauss/f?). Notice that these values of cAt 
are much smaller than the jet size [4] which justifies taking the magnetic field 
constant. 

Notice that according to Eqs. (45,49), the ratio a E /E u increases with 
energy. Thus, in the framework of the quasi-degenerate mass scheme, for very 
high energies and magnetic fields {i.e., for E v ^> 10 TeV and B ~ 10 9 Gauss), 
the spread of the wavepacket can overcome the wavepacket separation (see 
Eq. (16)). 

Our approximation of taking the momentum of the particle constant over 
At is valid only as long as AtE u 3> 1 . It is straightforward to check that for 

( ?. )( E » > 10 8 

l 10 7 Gauss 10 TeV T jet ' 

this condition is not fulfilled, so Eqs. (45,49) which were derived on the basis 
of this approximation are not valid. For 

' B :H^)(^)~10<, 



v 10 7 Gauss /y 10 TeV y T je t' 

in both muon and pion decay cases, oeIE v ~ 0.1 and the coherent correction 
to the spectrum can be significant. We will discuss this point in section 4.4. 



4 Applications 

4.1 Measuring the size of the wavepakets 

As we discussed in the previous sections the sizes of the neutrino wavepack- 
ets carry important information on conditions of neutrino production: on 
strength of the magnetic field, on radiation and matter density, on distance 
to the source. The sizes of the packets differ for neutrinos from the muon 
and pion decays. Can the sizes of the wavepackets of cosmic neutrinos be 
measured? Are they relevant for observations? 
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In Ref. [14], it has been shown that the density matrix of a "stationary" 
beam of particles is entirely determined by the energy spectrum and does 
not depend on the sizes of the wavepackets. As a result, all the "station- 
ary" beams composed of wavepackets of different sizes are equivalent and 
cannot be distinguished by measurement provided that they have the same 
energy spectra. However, in reality, even a macroscopically stationary beam 
composed of wavepackets of finite size is not stationary at the microscopic 
level. As pointed out in [14] , for a non- stationary beam, using time informa- 
tion in principle makes determining the size of the constituent wavepackets 
a possibility. 

Consider a burst of neutrinos emitted from a source at distance of L 
during time interval t b . If the emission is spherically symmetric, at the time 
of detection neutrinos will occupy a shell with a volume of AaiL 2 t h . The total 
number of emitted neutrinos is E burst /E u , where E burst is the total energy 
of the burst and E v is the average energy of neutrinos. The average time 
between two successive neutrinos passing through an area of S (perpendicular 
to the line of sight of the source) is 



[1 sec) 



V ^ J V Eburst 

L \ fl m 2 \ / E v \ /10 52 erg\ / t b 



100 MpcJ V S J \10 TeVy \ Eburst J V 1 sec / 
The time resolution of detector can be much better than At = 1 sec. As a 
result, the flux will appear as non- stationary at the microscopic scale. Sup- 
pose information on the arrival time of each wavepacket is provided. Then, 
by recording the detection time with precision better than At, one can in 
principle extract information on the size of wavepackets composing such a 
flux. If a detection takes place at time r from the arrival time of the "center" 
of a certain wavepacket, the wavepacket size has to be larger than cr. The 
question is whether there is any possibility to obtain the arrival time of each 
wavepacket. The answer is positive if, for example, the injection (production 
of pions) has some fine structure with 5t < cr x <^ At, associated to emission 
of the short 7— burst or gravitational waves; i.e., in this case, the arrival 
time of the "center" (or initial moment) of the neutrino wavepacket can be 
determined from the arrival time of 7's (or gravitational waves). Unfortu- 
nately, present models do not predict such strict time correlation (5t At) 
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between the 7 and v bursts. In more realistic situations when no information 
on spatial distribution of the wavepackets within the duration of the burst 
is available, even recording the detection time will not help us to determine 
the wavepacket size. 

The present situation argument is in line with the theorem proved in 
sec. D.3 of [13] which claims that for any flux composed of "short" wavepack- 
ets with a certain spatial distribution, it is possible to find an equivalent flux 
(i.e., having the same density matrix) composed of long wavepackets with 
width Az ^> 8z provided that the long wavepackets cover the same part of 
the space in which the short wavepackets are distributed. Although the sizes 
of the wavepackets cannot be practically determined, their effects can appear 
in the energy spectrum through the coherent broadening phenomenon that 
is discussed in sec. 4.4. 

4.2 Coherence of neutrinos from various cosmic sources 

In this section, we check if the condition of coherence loss can be satisfied 
for various sources of high energy neutrinos. The most powerful sources of 
neutrinos (Gamma Ray Burst sources and AGNs) are located at very far 
distances (L > 1 Mpc). From Eqs. (47,51), and for known values of Am 2 
we observe that even a relatively small magnetic field is enough to cause loss 
of coherence over such distances. The magnetic fields in the GRB sources 
[4] are of course larger than these bounds. Notice that the bound shown in 
Eq. (51) is independent of the energy while the bound in Eq. (47) becomes 
stronger with energy. Taking E v as high as 10 21 eV, from (47) we conclude 
that B larger than 0.1 milli-Gauss at the source is enough to destroy the 
coherence of the neutrinos from the pion decay over cosmological distances. 

Now let us discuss the cosmogenic neutrinos. Even though the baseline 
for these neutrinos is of cosmological scale, the coherence of the neutrinos 
produced by the muon decay is not lost because in the intergalactic area, 
where these neutrinos are produced, the magnetic field is too weak to satisfy 
Eq. (51). However, it is not out of question to have an intergalactic mag- 
netic field large enough to satisfy Eq. (47) for E v < EeV, and cause loss of 
coherence for the cosmogenic neutrino from the pion decay [16]. 

For neutrinos arriving from the center of galaxy (L ~10 kpc), to satisfy 
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Eq. (47) the magnetic field at the source has to be larger than 5 x 10~ 5 Gauss 
while to satisfy Eq. (51) the value of the magnetic field has to be larger than 
10 5 Gauss. If the magnetic field is between 5 x 1CT 5 Gauss and 10 5 Gauss, 
the part of the flux from the pion decay will lose its coherence but the other 
part from the muon decay will maintain its coherence. 

For "solar atmospheric neutrinos" (i.e., neutrinos produced via interac- 
tion of the cosmic ray particles with the atmosphere of the sun [18]), neither 
of the bounds in Eq. (47) and Eq. (51) is satisfied: Setting L = 5 x 1(T 6 pc 
(the distance between the Sun and the Earth), Tj et = 1, Am 2 = Am 2 atm and 
E v ~ 1 GeV, we find that in order to satisfy Eqs. (47) and (51) the magnetic 
field at the solar surface has to be larger than 10 5 Gauss and 10 11 Gauss, 
respectively; that is while, the magnetic field on the surface of the sun is less 
than 10 4 Gauss so the coherence is not destroyed. 

4.3 Coherence of neutrinos from terrestrial sources 

In this section, for comparison we discuss the coherence in terrestrial long 
baseline experiments with a strong magnetic field at the production site. 
Taking E u ~ 10 GeV, L ~ 2i? ffi , T jet = 1 and Am 2 = 2 x 1(T 3 eV 2 , from 
Eq. (47) we find that the magnetic field at the source has to be larger than 
10 9 Tesla to obtain complete loss of coherence. For Am 2 ~ 1 eV 2 or larger 
(as suggested by LSND; see [17] for recent status of such models), a magnetic 
field of 1000 Tesla at the source would be required to cause loss of coherence 
over a baseline of the Earth size. From Eq. (45), for E v ~ 30 GeV, Tj et = 1 
and B ~ 1 Tesla, we obtain <7e/E u ~ 10~ 6 . In practice, the energy resolution 
of the detector will be larger (worse) than this value. 

For a setup such as the muon storage ring of the neutrino factory the 
magnetic field is significant only at rounded edges of the rectangular track 
of the muons. The neutrinos reaching the detector are produced while the 
muon propagates through the side of the track that is aligned towards the 
detector. The length of the wavepacket of neutrino is therefore determined 
by the length of this side of the track, l str . Taking into account Lorentz 
contraction (see sec 3.1) we find that the length of the wavepacket emitted 
in the direction of the muon momentum is equal to 

&X — Istr-^-- (53) 
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Taking = 30 GeV and l s t r = 600 m, we obtain a x = 0.7 cm. From 
Eq. (14), we find that to break the coherence for the atmospheric mass split- 
ting the baseline has to be larger than ~ 0.1 Mpc. Obviously, the coherence 
of the neutrinos reaching the detector of the neutrino factory will be main- 
tained. 



4.4 Effects of coherent broadening 

Coherent broadening can lead to deformation of the energy spectrum. A 
well-known example is the coherent broadening of a monochromatic line in 
atomic spectroscopy. The effects smooth down any sharp structure of scale 
smaller than <je in the spectrum. If the coherent broadening is symmetric 
(i.e., if the broadening takes "a line" of energy E to (E — <7e/2, E + a E /2) ), 
the leading order correction to the spectrum (F) will be given by 



■ a 2 d 2 F 

F + — 

2 dE 2 



1 

N' 



where N is a normalization factor. However, if the broadening is asymmetric, 
i.e., if E is mapped onto (E — <7e(1/2 + a),E + o"e(1/2 — a)) with a ^ 0, the 
first order correction to the spectrum will be larger and linear in <je, i-e., 

^ dF 
F+ dE aEa 

In the conservative case of symmetric broadening, the effects on spectrum can 
be significant if oe ~ \(d 2 F / dE 2 ) / F]~ l l 2 , where F is the spectrum without 
taking into account the coherent effects. Notice that the effect takes place 
at production, so the distance between the source and detector does not 
play any role in this consideration. For cosmic neutrinos, the spectrum is 
predicted to follow a power law behavior: F oc E~ n with n ~ 2. For such a 
spectrum, the coherent broadening can be significant if oeIE v ~ 0.1. From 
Eqs. (33,34,35), we observe that the coherent broadening due to scattering 
is completely negligible. However, as discussed in section 3.3, for the decay 
in the magnetic field B ~ 10 9 GaussFj et (10 PeV/E„), the ratio a E /E can be 
of order of 0.1 (see Eqs. (45,49)). For lower values of the magnetic field and 
energy, the coherent broadening is irrelevant. As discussed in [6], in mildly 
relativistic jets with Tj et ~ 3 present in type Ib/c supernovae, the magnetic 
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field can be as high as 10 9 Gauss. In such an environment, the corrections 
to the spectrum due to the coherent broadening can be significant if E v > 
few 10 PeV. In principle, the energy loss through synchrotron radiation can 
hamper production of such high energy neutrinos. Notice, however that as 
shown in [26] through the same mechanism that protons are accelerated in 
the magnetic field, muons and charged mesons can also be accelerated. As a 
result, reaching energies above PeV is possible [26]. Since the dependence of 
<j e on energy is non-linear (see Eqs. (45,49)), the coherent effects can change 
the shape of the energy spectrum. 

4.5 Effects of incomplete averaging 

As we discussed in sec. 4.2, in some special cases, the coherence of cos- 
mic neutrinos will be maintained. In case that coherence is not lost, the 
probability of neutrino oscillation can be represented as 

P = P + P mt , (54) 

where P is the average probability and Pj ni is the interference term which in 
the two neutrino scheme can be written as 

P int = D cos <I>(E,L) . (55) 

P and D are independent of the energy of the neutrinos, E, and the baseline, 
L. The phase of oscillations in vacuum is defined in Eq. (5). In the case of 
three neutrino oscillation, there will be three interference terms associated 
to three different Am 2 . 

The oscillatory terms are in practice averaged to zero for the following 
three reasons: (i) the integration of signal over the energy because of the 
finite energy resolution of the detector; (ii) the integration over the region of 
the production; and (iii) the summation over different cosmic sources. The 
third reason is specially inevitable. For example, from a single gamma burst 
source at z ~ 1, we expect only 10 _1 — 10 events at a 1 km 3 scale detector 
[4]. Fortunately, about ~1000 GRB are observed every year which can yield 
at least around a few hundred neutrino events per year at ICECUBE. When 
we combine the data from different GRB sources, averaging over baselines 
takes place. 
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For typical values of the cosmic neutrino parameters, we obtain the phase 



,„ / Am 2 \ / L \ /10 TeV\ , . 

Considering the fact that the phases are very large, <^ > 1, we expect the 
oscillatory term to average to zero. In the following, we discuss the error 
from neglecting these terms. We discuss the procedure of averaging for two 
neutrino mixing where only one interference term exists. Generalization of 
our results to the three neutrino case is straightforward. For simplicity and 
definiteness, we will consider averaging due to integration over the neutrino 
energy. Integration over source volume and summation over different sources 
do not produce new features and can be considered in a similar way. The 
"average" number of events of different types in the detector can be written 
as 

(57) 

where E th is the detection threshold and 

dN dN°{E) 



dE ' dE 



-P(E,L). 



Here dN°(E)/dE gives the average number of events with neutrino energy 
E which would be observed if there was no oscillation [i.e., if P{y a — > v$) = 

Let us divide the whole detectable energy range into very small intervals 
AE(<ti AEl) over which the dependence of dN/dE on energy can be ne- 
glected. The contribution of each interval to the total number is given by 
a Poisson distribution whose average and variance are both determined by 
(dN/dE)AE. According to the "central limit theorem" [19], the distribution 
of the total number is given by a Gaussian whose mean and variance (<r 2 ) 
are both given by a sum (an integral) over (dN/dE): 

{(N - <iV» 2 > = (N) , 

where (N) is given by Eq. (57). Because of the uncertainty in the mea- 
surement of L, we cannot in practice determine the value of cos<p(E, L) and 
therefore that of dN/dE. However, regardless of the value of <p(E,L), the 
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integral of the oscillatory terms (see Eq. (57)) vanishes up to a small thresh- 
old correction: 0[AEl(cIN / dE)\ t h] <C VN. Thus, in the analysis of cosmic 
neutrinos, we can safely drop the oscillatory terms. 

As shown in [13], from the observational point of view, the effect of aver- 
aging of oscillation is equivalent to the effect of coherence loss when signals 
from separated wavepackets sum up in the detector. As discussed, the vari- 
ance of the distribution of the number of observed events is also the same for 
both cases and is determined by the root of the number of events regardless 
of the shape of the neutrino spectrum. 

4.6 Kinematical decoherence versus quantum decoher- 
ence 

Throughout this paper, we have assumed that the evolution of the neutrino 
states from the source to the Earth is given by the standard quantum mechan- 
ics. As is well-known, within the quantum mechanics, a pure state remains 
pure forever. Quantum gravity is expected to modify the evolution of the 
quantum states [20, 21]. As long as the modifications preserve unitarity and 
the average energy and moreover respect the second law of thermodynamics, 
the damping factors appear only in front of the oscillatory terms [22, 23]: 

P(u a -> up) ~ P{y a -> up) ~ \U al \ 2 \U pi \ 2 + \U a2 \ 2 \Up 2 \ 2 + \U a3 \ 2 \Up 3 \ 2 

+2Re[U; i U al Up 2 U* a2 ] cos A 12 e~ L ™ + 2Re[U*p 1 U al Up 3 U* a3 ] cos A 13 e- L ^ 
+2Re[^ 2 f/ Q2 ^3^3]cosA 32 e- i ^ , 

where Ay = Am- L/(2E V ). Thus, in the case of the cosmic neutrinos, for 
which either the interference terms disappear or their effects average to zero, 
the new damping factors are irrelevant. To have a deviation from the stan- 
dard picture some of the assumptions that went into this conclusion (e.g., 
conservation of the average energy) have to be relaxed [24]. Our results are 
in accord with [25] . 

5 Conclusions and Summary 

In this paper, we have explored various aspects of the coherence of the cosmic 
neutrinos with special emphasis on the neutrino flux in the (1-100) TeV 
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energy range. Since the three neutrino mass eigenstates composing a single 
flavor state have different velocities, they spatially separate on their way 
to the detectors (see Eqs. (8,12)) provided that the source is far enough. 
Their effects in a detector cease to interfere and the oscillatory behavior of 
the probabilities disappears. The degree of separation depends on the size 
of the wavepacket, which in turn, is determined by the properties of the 
environment at the source. 

We have evaluated the length of the neutrino wavepackets produced by 
the pion and muon in various environments. In an interaction-free medium, 
the wavepacket sizes in the observer frame from both the pion and muon 
decay are determined by the lifetime of the parent particle, the energy of the 
parent particle and direction of neutrino with respect to the momentum of 
the decaying particle. The size of the wavepacket does not depend on the 
feature of the decay (two-body vs. three-body decay). In contrast, when 
the parent particles undergo interactions before decay, the wavepacket sizes 
become dramatically sensitive to the shape of spectrum. That is despite the 
fact that the electromagnetic interactions of the muon and pion are practi- 
cally the same, the neutrino wavepacket sizes and their dependence on the 
parameters become completely different for the pion and muon cases. 

We have found that in an interaction-free environment, the wavepacket 
length is large enough to maintain coherence even for neutrinos traveling 
over cosmological distances. So the high energy cosmogenic neutrinos should 
arrive in coherent states 

Scattering of the parent charged particle off the particles present in the 
medium decreases the length of the produced neutrino wavepackets. We 
have evaluated the amount of shortening of the wavepacket in a radiation 
dominated environment such as the one described in [4] as the GRB source. 
We have found that in the case of neutrinos produced by the pion, scatter- 
ing can cause complete loss of coherence (even for Am 2 sol ) over cosmological 
distances. However, in the case of neutrinos produced by the muon with 
E v ~ 100 TeV, the effects of scattering are too small to cause coherence loss. 
In contrast, we have shown that even a tiny magnetic field at the production 
region can considerably widen the width of the neutrino wavepacket. A mag- 
netic field as low as a few ten Gauss at the source is enough for effective loss 
of coherence (see Eqs. (47,51)). Based on the above consideration, we have 
concluded that because of the sizeable magnetic fields, all the neutrinos from 
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remote point sources such as GRBs and AGNs lose their coherence before 
reaching the Earth. We have found that the dependence of the wavepacket 
size on the magnetic field and energy is quite different for neutrinos from 
the pion and muon decay. As a result, it is conceivable to have a situation 
that neutrinos from pion decay lose their coherence but those from the muon 
decay maintain theirs. Such a situation can be realized for cosmogenic neu- 
trinos as well as for "solar atmospheric neutrinos" . For all realistic terrestrial 
setups the coherence is maintained. 

We have found that the coherent broadening can be important for 

Gauss r jet (10 PeV/E u ). 

Although such condition seems quite extreme, its realization is not out of 
question [6, 26]. The dependence of oe on energy is non-linear [see Eqs. (45,49)] 
As a result, the coherent effects can change the shape of the energy spectrum 
of cosmic neutrinos. 

We have studied the case in which the coherence is maintained on the 
way from the source to detector. We have found that the variance of the 
distribution of the number of events is the same as the variance that we 
would expect in the absence of the oscillatory terms. As a result, neglecting 
the oscillatory terms does not induce a new source of error in determination 
of the neutrino parameters. 

We have also discussed the effects of a deviation of the evolution of neu- 
trino states from that in the standard quantum mechanics as one would ex- 
pect in the framework of the quantum gravity. We have shown that as long 
as the deviation preserves unitarity and the average energy and moreover 
respects the second law of thermodynamics, the neutrino flavor composition 
at the detector is the same as the one without the deviation. To obtain a 
different composition, some of these assumptions have to be relaxed. 

The neutrino pattern at the detector turns out to be the same for all 
three cases that we have studied in the present paper (i.e., maintenance and 
loss of coherence and a deviation from standard quantum mechanics) and 
as a result, these three cases cannot be distinguished by observation. Other 
possible applications of the results will be considered elsewhere. 
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Appendix 



In the following, we study the condition for the coherence loss in the energy- 
momentum space. Let us consider a neutrino field of flavor a characterized 
by the following wavepacket 

\v a ) = J f(p)\a;p)d 3 p, 

where f(p) is a function of the momentum which is nonzero only in a small 
interval with width a p . The amplitude of finding the neutrino with flavor f3, 
up, at point x and at time t is 




where p = \p\. Let us denote the spatial resolution by AX D . The probability 
of finding vp in a volume of V = (AX D ) 3 at time t = L (where L is the 
distance between the source and the detector) is 

P(v a -> vp) = / \(up;x,L\u a )\ 2 d 3 x . (58) 
Jv 

For any realistic detector, AX D is far larger than the size of the wavepacket. 
As a result, we can replace the integration over x with a delta function: 

f(p) is a relatively smooth function over an interval of size a p . For (m| — 
mf) satisfying the condition of loss of coherence (8), the integration of the 
interference terms over p vanishes and we therefore obtain 

P(^ a ^^)=Y t \U fii \ 2 \U 0li \ 2 . 

i 

Notice that this is the probability for a single neutrino rather than the average 
probability for several neutrinos states. 

As discussed in [13], even if the condition (8) is fulfilled, for a hypothet- 
ical detector with energy resolution, AEp, smaller than a p , the interference 
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terms can become relevant provided that (Am 2 L/p)(AEf)/p) ~ 1. Let us 
discuss how the difference between the two cases can be understood within 
the formalism discussed in this appendix. To calculate the total probability 
of detecting up at the detector, we first calculated the probability of appear- 
ance of U/3 at a certain point (x) of the space and then integrated over x. 
The rationale behind this treatment was that in a realistic neutrino detec- 
tor, the neutrinos are detected through the electroweak interactions whose 
range (i.e., ~ 1CT 16 cm) is much smaller than the typical size of the 
wavepacket. In terms of quantum mechanics, to a good approximation, the 
eigenstates of a realistic detection operator are \vp;x,t). In contrast, the 
eigenstates of a hypothetical detector with an energy resolution of AEp are 
states described by \vp,D) = f gD{p)\vp] p)d 3 p \v$\ x, L)), where goip) 
is a function of momentum with width AE D . Thus, with such a detector, 
the total probability is given by an integral over \ {vp\ -D|z/ a )| 2 rather than by 
Eq. (58). 
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